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Abstract 

The critical properties of the Ziff-Gulari-Barshad (ZGB) model with the addition of long-range 
reactivity [G. H. Chan and P. A. Rikvold, Phys. Rev. E 91 , 012103 (2015)] are further investigated. 
The scaling behaviors of the order parameter, susceptibility, and correlation length provide addi¬ 
tional evidence that the universality class of the ZGB system changes from the two-dimensional 
Ising class to the mean-field class with the addition of even a weak long-range reactivity mechanism. 


PACS numbers: 05.50.-tq,64.60.Ht,82.65.-tr,82.20.Wt 


*Electronic address: prikvold@fsu.edu 


1 



I. MODEL AND MOTIVATION 


In 1986, Ziff, Gulari, and Barshad introdnced a lattice-gas model (known as the ZGB 
model) to simnlate the noneqnilibrium process of the oxidation of carbon monoxide on a 
Pt catalyst j^. The catalyst snrface is modeled as a L x L sqnare lattice, and a random 
site is chosen times in each time step (MGSS, Monte Garlo step per site). Adsorption 
of the gas species can take place only if the chosen site is empty. GO is adsorbed with 
probability y, which is proportional to its partial pressnre in the supplied gas. Adsorption 
of oxygen is attempted with probability 1 — y, but is successful only if a randomly chosen 
nearest-neighbor site is also empty. O 2 then dissociates into two adsorbed O atoms. After 
the adsorption, the adsorbed species checks its four nearest-neighbor sites. If the opposite 
species is found, the pair forms carbon dioxide (GO 2 ) and leaves the surface immediately. It 
has further been noticed in experiments that adsorbed species can desorb from the surface 
without reacting, and that the desorption rate of carbon monoxide is much higher than that 
of oxygen j^. Therefore, the desorption rate of GO (k) is often added to the model as a 
second control parameter, in addition to the GO partial pressure (y) [l, 6 |. In this work, we 
choose desorption with probability k. The desorption rate can be considered as a proxy for 
temperature as thermal fluctuations would make the adsorbed particle break its bond with 
the catalyst surface more easily. 

If the GO desorption rate {k) is not very high, simulations show that the catalyst surface 
can be in one of three phases, depending on the value of the GO partial pressure (y): 1. GO 
poisoned phase, which has GO covering most of the surface with some scattered empty sites. 
It is observed when y is large. 2. 0-poisoned phase, which has oxygen covering the whole 
surface when the GO partial pressure [y) is sufficiently low (not observed in experiments). 
3. mixed phase, in which the catalyst surface is covered by a mixture of O, GO, and empty 
sites. It is observed for intermediate values of y. The order parameter of the system is the 
GO coverage ( 6 ^ 00 ), which is dehned as the proportion of the lattice sites occupied by GO. 
A hrst-order phase transition line separates the mixed phase and the GO poisoned phase. 
Similar to an equilibrium lattice-gas system, moving along this hrst-order transition line in 
the phase diagram eventually leads to a critical point. It has been found that this critical 
point belongs to the two-dimensional equilibrium Ising universality class 6 |. 


Recent studies of adding non-zero long-range interactions in an equilibrium Ising ferro- 
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magnet showed that the universality class of the critical point changed abruptly from 
Ising to mean-held. Motivated by this result, we recently modihed the ZGB model by in¬ 
troducing an adjustable probability (a) that an O atom and a CO molecule adsorbed far 
apart on the surface can react to form CO 2 and desorb [^. This move is attempted only if 
the adsorbed species cannot hnd an opposite species within its four nearest-neighbor sites. 
If the two sites contain opposite species, CO 2 will be formed and the two sites will become 
empty. Details of the implementation are shown in Fig. [H 

While our recent work has shown that the universality class of the critical point 
changes from the Ising class to the mean-held class using fourth-order cumulants, here we 
look at the huctuations of the order parameter (CO coverage) at the critical point through 
two quantities, the ’susceptibility’ and the time average of the absolute diherence of the CO 
coverage from its mean value. The results further conhrm the change of the universality 
class of the critical point from Ising to mean-held upon introducing long-range reactivity. 


II. RESULTS 

Figure [2](a) shows a Ln-Ln scaling relation between the time average of the absolute 
diherence of the CO coverage from its mean ((|6*co,c — (6*co,c)|)) at the critical point, and 
the system size (L). The critical points were found in our recent work a. cu^ulant 

crossing. The negatives of the slopes of these lines represent the critical exponent ratio (ijv 
and are plotted in Fig. [2](b). Without long-range reactivity (a = 0), we obtain 0.0977 ± 
0.0007, which is close to the exact Ising value of /5 = 1/8 and v = 1. For the case with 
long-range reactivity (a > 0), the weighted average is 0.564 ± 0.008, consistent with the 
exact mean-held result of = 1/2 and v = 1. 

We dehne a nonequilibrium analog of equilibrium magnetic susceptibility or huid com¬ 
pressibility as 

XL = L\{elo,L) - {OcO.L?) (1) 

for a. L X L system, which for simplicity we call ’susceptibility.’ It is well suited to measure 
the strong huctuations of the order parameter near the nonequilibrium critical point jd]. 
Figure. [3](a) shows the natural logarithm of the susceptibility {xl) vs the natural logarithm 
of the system size (L). The slopes of the htting lines represent the critical exponent ratio 7 /z/ 
and are shown in Fig. [3](b). For a = 0, we obtain 7 /z/ = 1.7967 ± 0.0001, which is close 
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to the exact Ising value of 7/4. For a > 0, we get a weighted average of 7 / 1 / = 0.92 ± 0.01, 
which is close to the exact mean-held value of 1 . 

We dehne the CO disconnected correlation function as c(r) = (criaj), where Uj is 1 if site 
i is occupied by CO and is 0 otherwise, r is the distance between site i and site j, and the 
spatial average is taken along the horizontal and vertical directions. The critical correlation 
length was estimated by integration as 


- (c(L/2)}]rdr 




- (c(L/2)}]dr 


and j^. 


( 2 ) 


as m [ 0 |j ana [^j. In the study of an Ising model with long-range interactions of strength a 
j^, Nakada et al. suggested the scaling relation ^{L) = They plotted ^{L)/L 

against using the Ising value of 1^/7 = 4/7 and verihed that the data points all fell onto 


one curve When a was not too small, the log-log plot of the graph had a slope of —1. 
Here we tried the same procedure with our long-range reactivity a, also using 1^/7 = 4/7, 
and obtained Fig. ID^a). Although the data points do not he perfectly on one curve, in a 
Ln-Ln plot, we also obtain a straight line with slope close to —1, as shown in Fig. IH^b). 
Using zz /7 = 1/1.7967 as obtained in Fig. IH gives nearly exactly the same graphs, just with 
the slope in (b) changed from — 1.02 to —1.03. So it is justihed to claim 1^/7 = 4/7. 


B. 


In our recent work l2|, we also tried to obtain 1^/7 by looking at the scaling relation 
between the change in critical desorption rate with the change in long-range reactivity, as 
well as the scaling relation between change in critical CO partial pressure with the change in 
long-range reactivity. We obtained fcc,oo(o) — A;c,oo(0) oc and 2 /c, 00 (a) — i/c,oo(0) oc 
and expected that one of these values should be 1^/7 of the a = 0 case, as we know that 
for the 2D equilibrium Ising model with long-range interaction of strength a, the change 
in critical temperature against the change in long-range interaction strength a is known to 
behave as Tc { a , 00) — Tc(0, 00) oc with 1^/7 obtained from the a = 0 case, i.e. 4/7 j^. 

The index values of 0.499 and 0.448, both deviate signihcantly from 4/7. We therefore 
suggested in that a might not be linearly related to a, but rather as a oc a^'^. However, 
Fig. mb) suggests that the relation for the equilibrium system, ^{L) = still holds 

for our model with a replaced by a, which means a oc a. Both Fig. [3](b) and Fig. lU^b) suggest 
that Z//7 = 4/7 for a = 0. Thus it seems that the indices obtained from — fcc,oo(0) 

vs a and i/c, 00 (a) — i/c,oo(0) vs a are not Z//7. We currently have no satisfactory explanation 
for this result. 
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III. SUMMARY 


We present numerical evidence that the addition of non-zero long-range reactivity to the 
ZGB model changes the universality class of the critical point from Ising to mean-held. 
Estimates for the exponent ratios (ijv and 7 / 1 / are obtained and compared to the value for 
7 / 1 / obtained in our previous published work j^. 

This work was supported in part by NSF grant No. DMR-1104829. 
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FIG. 1: Flow chart for the reaction process. The algorithm is based on 
ZGB model with CO desorption. The framed region contains the added 
strength a. (From [3]). 
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FIG. 2: (a) Natural logarithm of the time average of the difference of the CO coverage from its 

mean, shown vs the natural logarithm of the size of the system for L = 60, 100, 160, and 240. 10^ 
MCSS were used for non-zero long-range reactivity (a = 0.1, 0.3,0.5, 0.7 and 1), and 10® MCSS 
were used without long-range reactivity (a = 0), when the system was near the corresponding 
critical point (e.g, the data point for L = 100 was taken at the 100/60 cumulant crossing point). 
The negative of the slopes of these lines are the critical exponent ratio fijv and are plotted in (b). 
For a = 0, it gives 0.0977±0.0007, which is close to the exact Ising value of 1/8 = 0.125. For a > 0, 
the weighted average value is 0.564 ± 0.008 (dashed line), which is close to the exact mean-field 
value of 1/2 (dotted line). 


7 










FIG. 3: (a) Natural logarithm of the susceptibility (xl), shown vs the natural logarithm of the 

system size (L) for L = 60, 100, 160, and 240 near their corresponding critical points with o = 0, 
0.1, 0.3, 0.5, 0.7, and 1 . 5 x 10^ MCSS were used for a > 0 and 10® MCSS for a = 0. The slopes of 
the fitting lines for different long-range reactivities (a) are the critical exponent ratio plotted 
in (b). For the Ising case (a = 0), we obtain 7 / 1 / = 1.7967 ± 0.0001, which is close to the exact 
Ising value of 7/4. For the mean-field case (a > 0), we get a weighted average of 7 / 1 / = 0.92 ± 0.01 
(dashed line), which is close to the exact mean-field value of 1 (dotted line). 



FIG. 4: (a) Scaling plot of the correlation length at the corresponding critical points with a 

scaling index of 4/7 for a > 0, using lO’^ MCSS. Scaling index values of 0.435,0.448,0.499,0.45 and 
(4/7)/1.3 were also tried but not shown, as 4/7 gives the best result, i.e., the data points fall most 
closely onto one curve. A log-log plot of the results in (a) are shown in (b). The data points he 
along a straight line of slope close to —1, which justifies the use of 4/7 as the scaling index. 
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